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Abstract 

The parabolic approximation is developed for high energy charged particles scatter- 
ing in a bent crystal with variable curvature. The general form of parabolic equation 
is received for atomic chains located along coordinate axis of orthogonal curvilinear 
coordinate system. 
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1 Introduction 



The parabolic equation was introduced by Leontovich and Fock [1] to de- 
scribe the scattering of the wave in a cone centered on a direction of wave 
propagation. Lervig, Lindhard and Nielsen [2] tried to use the similar time de- 
pendent Schrodinger equation for quantum treatment the directional effects 
of energetic charged particles in the crystals. The work [2] indicated that an 
asymptotic expansion of the wave equation with large parameter p * R 3> 1 (p 
- the incident momentum of particle, R - the atomic screen radius) should not 
be done in the direction of p, but in the direction of crystal axes (planes), which 
is slightly different from the direction of the particle. The work [3] shown that 
accurate calculations of potential scattering in the field of a straight atomic 
chain can be done as a double asymptotic expansion in terms of two large pa- 
rameters of the scattering: the length of the atomic chain L x /R 3> 1 and the 
component of the incident momentum of particle along the chain p x * R^> 1. 
It is important to note [3,4] that this double expansion gives the parabolic 
equation as well as the model of continuum potential which is central idea 
of channeling effect. The expansion along the axis of symmetry is the major 
difference between original parabolic equation [1] and its applications for the 
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scattering in the crystals. The parabolic equation method is also more general 
than sudden collision approximation [5] or so called eikonal-type approxima- 
tions [6,7], and semi-classical corrections to the eikonal scattering amplitude 
can be derived from it. 

In this report we receive the parabolic equations using the symmetry of the 
crystals described by orthogonal curvilinear coordinate system. On local level, 
such systems looks as Cartesian coordinate systems. However, there is a differ- 
ence between the standard parabolic equation [1] and the parabolic equations 
in curvilinear coordinate systems considered here. 



2 Parabolic equations in cylindrical symmetry 



The parabolic approximation allows to construct a family of solutions to the 
Schrodinger wave equation which are close to the plane wave along some di- 
rection. If the potential of scattering U has a spherical symmetry, the only 
possible direction of expansion is the direction of incidence p of the particle. 
Denoting this direction as x, the wave function \1> satisfying to the stationary 
Schrodinger equation 

(A+p 2 -V)^ = 0, (1) 
and can be presented as 



* = exp(ipx)®(r), (2) 

where V = 2MU, M is the mass of particle (H — c — 1), <&(r is very slowly 
varying in the direction x function in comparison to the expiipx). Substitution 
of Eq. (2) into (1) and neglecting of the term d 2 $/dx 2 yields a well known 
parabolic equation [1,8] 



2ip d^/dx + (A ± - V)<$> = 0, (3) 

with transverse Laplacian Aj_ = d 2 /dy 2 + d 2 /dz 2 . When the scattering is 
studied on potential centers which constitute a straight atomic chain (or 
plane), and the momentum of particle is not parallel to the axis x of chain, 
pr = p x x + p±f± , the the wave function \I/ should reflect the translation 
symmetry along the direction of the chain 



\& = exp{ip x x)<&{f). 



(4) 



Fig. 1. Atomic chain along (a)the circumference and (b) radius of circle 
Then the parabolic equation becomes 

2ip x d$/dx + (A ± + p 2 - p 2 x - V)$ = 0. (5) 

Now we turn to the parabolic equation for bent atomic chain. First, consider a 
circlular bent with a constant curvature k = 1/Ri and parametric presentation 
of a chain along this circle x = Ri sm((p) , y = i?;(cos(0) — l),z = z (see 
Fig. 1(a)). If we take a cylindrical coordinate system (p,<j),z), then there are 
the conjugate momentums of free particle (pp,P4>,Pz), such that 



pr = p pP + P4>( f) + p z z. (6) 

Based on this presentation of phase, we can derive three different parabolic 
equations, each one for the propagation of the particle along one curvilinear 
coordinate when two others are constants. Parabolic equation along the axis 
z lookes the same as 5, so we consider two cases in Fig. 1. For the motion 
Fig. 1(a) the wave function can be in the form 



V = exp(ip <j> <i>)$(p,<l),z). (7) 
Using Laplacian in cylindrical coordinates 



A —— — -— — 

p 2 d<f) 2 dp 2 p dp dz 2 ' 



we obtaine the eq. 



with transverse 2D Laplacian Aj_ = d 2 /dp 2 + 1/pd/dp + d 2 /dz 2 . The next 
step is to neglect the small term 1/ p 2 d 2 §/d<f) 2 in range of high energies, and 
we have the final result 
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W* + ^ +p2 -7- v) ^°- (10) 

In the second case (Fig. 1(b)), the wave function can be written as 

* = exp(ip p p)$(p, 0, z), (11) 

and doing the similar procedure, one can get the parabolic equation for the 
radial motion 



KPPJT + -^ + (^±+P 2 -P 2 P + ^-V)^ = 0. (12) 
op p dp P 

However, using the assumption p p p ^> 1 the terms 1 / p<9$ / dp and ip p / p can 
be disregarded in Eq.(12) and parabolic equation becomes the same as Eq.(5) 
for the straight chain. 



3 Parabolic equations in orthogonal curvilinear coordinate system 



Now we derive the general form the parabolic equation. Assume an atomic 
chain is located along the curvilinear axis q± of some orthogonal curvilinear 
coordinate system (qi , q 2 , q 3 ). The Laplacian in this coordinate system can be 
written in the form 



1 ( d ,h 2 h 3 d s d d s , d t h x h 2 d ,\ 

{IS) 



hih 2 h 3 \dq! hi dqi dq 2 h 2 dq 2 dq 3 h 3 dq 3 J 

where hi = ^Jgu is Lame coefficients; gij is metric tensor of the space, gij = 
for i ^ j. For example, in the cylindrical coordinate system considered in 
Sec. 2 hi = l,h 2 = p,h 3 = 1 and , in the elliptic coordinate system x = 
Ri cos(g 2 ) cosh(g!), y = R[ sin(g 2 ) cosh(gx), z = q 3 Lame coefficients are hi = 
h 2 = R[^J cosh(qi) 2 — cos(qi) 2 , h 3 = 1 . The wave function may be presented 
as (obviously, due to symmetry there is no difference what coordinate qi can 
be chosen) 



^ = exp(ip qi qi)$(qi, q 2 , q 3 ). (14) 

The possibility of such presentation is based on relation between direct (91, q 2 , q 3 ) 
and duai(p qi ,p q2 ,p q:j ) spaces 



pr = p qi qi + P Q2 q2 + P qi qz- (15) 
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Acting in the same manner as in Sec. 2 , we can derive three different parabolic 
equations, each one for propagation of the particle along one curvilinear co- 
ordinate when two others are constants. We proceed with the wave function 
(14). Substituting (14) in wave equation (1) with Laplacian (13), the accurate 
equation for reduced wave function $ can be written as 



2i Pqi <9$ + 1 d h 2 h 3 d<$> + 1_#^> + 
hi dqi hih 2 h 3 dqi hi dqi h\ dqj 

where the tranverse Laplacian is 



A - = ^rCf <nrir > + ^ < 17 > 

hih 2 h 3 \dq 2 h 2 dq 2 dq 3 h 3 dq 3 J 

Now since the second derivative p-^r in (16) has relatively small value, the 
parabolic equation in orthogonal curvilinear coordinates becomes 



2tp^d^_ + ^ d h 2 h 3 d$ | 

hf dqi hih 2 h 3 dqi hi dqi 

+(A ± +p 2 - % + ^^^(^) " V)* = 0. (18) 
hf hih 2 h 3 oq 1 hi 

This is the most general form of parabolic equation, and all results received 
in Sec. 2 can easily be deduced from it. Indeed, in Cartesian coordinates qi = 
x, q 2 = y, q 3 = z and h\ — h 2 — h 3 — 1, now from (18) it follows the Eq.(5). In 
cylindrical coordinates qi = p,q 2 = (fi,q 3 = z and hi — 1, h 2 — p, h 3 — 1. From 
(18) it follows the Eq.(12). Parabolic equations along the other axes q 2) q 3 look 
exactly the same as (18) with cyclic substitution of indices. For example, the 
parabolic equation along the axis q 2 will have the form 



2t Pq2 d$ + 1 d hih 3 8$ + 
h\ dq 2 hih 2 h 3 dq 2 h 2 dq 2 

+(A ± +p 2 - § + ~ V)* = 0. (19) 

h 2 hih 2 h 3 dq 2 h 2 

Now it is straightforward to see that in cylindrical coordinates Eq.(19) has the 
form of the Eq.(10). 
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